Priedel-like Oscillations from Interstitial Iron in Superconducting Fei+yTeo.62Seo.38 
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Using polarized and unpolarized neutron scattering we show that interstitial Fe in supercon- 
ducting Fei+yTei-iSe^ induces a magnetic Friedel-like oscillation that diflracts at Qj^ — (|0) and 
involves >50 neighboring Fe sites. The interstitial > 2 fiB moment is surrounded by compensating 
ferromagnetic four spin clusters that may seed double stripe ordering in Fei+j,Te. A semi-metallic 
5-band model with (^|) Fermi surface nesting and four fold symmetric super-exchange between 
interstitial Fe and two in-plane nearest neighbors largely accounts for the observed diffraction. 



While superconducting Fei+j^Tei„2,Sea; shares band 
structure, Fermi surface m, and a spin resonance with 
Fe pnictide superconductors 0-11], the parent magnetic 
structures are surprisingly different. Fig.[lja) depicts the 
distinct magnetic unit cells with single striped order for 
122 arsenides (q,„ = i^jk)) Hi versus double stripes for 
Fei+j,Te (q„ = (^,0)) fff,^. In this letter we show that 
short range ordered glassy magnetism at (^,0) in super- 
conducting Fei+yTei-xSex {x — 0.38) arises from mag- 
netic Friedel-like oscillations surrounding interstitial Fe. 
A critical role of interstitial iron to stabilize the lamellar 
structure Q , enhance magnetism 10| , and reduce the su- 
perconducting volume fraction [11 1 was previously noted. 
Our results provide a quantitative microscopic view of the 
pivotal magnetic polaron. 

We used three co-aligned Fei+yTeo.62Seo.38 single crys- 
tals with total mass sa 20 g and y = 0.01(2) determined 
by EDX. Grown by a flux method [TT'|, the samples are 
tetragonal (space group P4/nmm) with low temperature 
(T) lattice parameters a = 3.791 A and c = 6.023 A. 
Magnetization and specific heat measurements yielded 
Tc = 14.0(2) K and a superconducting volume fraction 
of 92.9(7)% and 83(1)% respectively (Fig. [Sjb)). 

Neutron scattering was performed using the Multi Axis 
Crystal Spectrometer (MACS) at NCNR, where a 20 
MW reactor, liquid H2 moderator, and doubly focussing 
PG(002) monochromator provided an incident beam flux 
of 2.0 X 10^ n/cm^/s [12]. Twenty detection channels 
permitted mapping of elastic scattering throughout the 
(HKO) and (HOL) reciprocal lattice planes with Ei — 
Ef ^ 3.6 meV (Fig. ^a.)-{h)). High T measurements 
(T = 25 K) provided background to cancel the domi- 
nant elastic nuclear scattering. Polarized neutrons were 
used to establish the magnetic origin and polarization of 
the scattering. Spin polarized '^He gas held in glass cells 
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FIG. 1. (a) Fe-plane magnetic order in the 122 and 11 parent 
compounds, (b) Half unit cell of Fei+yTei-^Sea; showing the 
location of interstitial Fe in orange (Fei). 



within a vertical solenoid concentric with the sample ro- 
tation axis was used to select the vertical component of 
neutron spin before and after detected scattering events 
13]. The 5 meV flipping ratio was typically 56 and 8.4 
for Bragg scattering from AI2O3 and Fei+yTeo.62Seo.38 
respectively. The corresponding sample depolarization 
factor of 0.825 was T— independent between 4 K and 30 
K. A channel mixing and transmission correction for time 
dependent '^He polarization (r w 60 — 90 hours) - averag- 
ing 60 (42) for the non-spin-flip (spin-flip) channel - was 
applied to T— difference data. The measurement pro- 
tocol ensured less than 5% effect of varying cell trans- 
mission on T— difference data. Absolute normalization 
of the unpolarized scattering cross section was obtained 
through comparison to acoustic phonon scattering and 
checked against incoherent elastic scattering from vana- 
dium. The polarized beam conflguration was calibrated 
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FIG. 2. (a) Constant hu — slice showing the difTerence be- 
tween neutron scattering intensity at T = 1.6 K and T = 25 K 
in the (HKO) and (b) (HOL) scattering planes. Inversion 
symmetry was imposed. Features near the origin, (0,0, ±1) 
Bragg peaks, and around the perimeter in (b) arise from in- 
tense nuclear scattering, (c)-(d) Calculated intensity distribu- 
tion for a 4-foId symmetric spin cluster surrounding intersti- 
tial Fe. (e)-(f) Calculated intensity for an interstitial Fe-site 
exchange coupled to a five d-orbital model. 



to the unpolarized configuration through incoherent elas- 
tic scattering from the sample. 

Fig. [Ha)-(b) show the wave vector dependence of the 
difference between elastic scattering at T = 1.6 K and 
T = 25 K. The rod like nature of scattering in the (HOL) 
plane (Fig. ^b)) indicates quasi-2D correlations. Ne- 
glecting the interstitial site, Fei+yTei-ajSe^; has only one 
Fe site per primitive unit cell. The wave vector depen- 
dence of magnetic neutron scattering associated with the 
periodic structure therefore must repeat in each Brillouin 
zone. This implies - modulo the magnetic form factor and 
polarization factor - that the intensity for L = ±| should 
match that a,t L — A possible explanation for the 
reduced intensity at L = ±| (Fig. [l{b)) is an uniaxial 
spin configuration ||c, which would imply ASz = in the 
corresponding neutron scattering cross section. 

Fig. [3] shows the Sz resolved T— difference scattering 
versus Q = (0.535, k, 0). The magnitude of the polarized 
scattering cross section is consistent with the unpolarized 
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FIG. 3. Polarized neutron scattering measured along 
Q=(0.535,K,0). AJc = [/(4 K)-/(25 K)]^ is the T-difference 
intensity without energy analysis following ^He cell transmis- 
sion correction. Red open diamonds are non-spin-flip data 
and blue open squares are spin flip data. The 1.4 mT guide 
field was perpendicular to Q and parallel to c. 



configuration (Fig. [5]). The peak in the spin- flip channel 
proves that part of the scattering cross section is mag- 
netic. Assuming the non-spin-flip T— difference intensity 
is also magnetic, the intensity ratio of 0.67(12) between 
the ASz = 1 and ASz = channels, implies that same 
ratio between the in- and out-of-plane components of the 
spin correlation function. This ratio is too large for the 
corresponding polarization factor to account for the re- 
duced intensity at L = ±|. Likewise, in the (HKO) plane 
the elastic magnetic scattering, which comprises four tri- 
angle shaped features at (±i,0) and (0,±i) is strongly 
suppressed in the adjoining Brillouin zones (Fig. [2][a)). 

Because the polarization and form factors for magnetic 
neutron scattering cannot account for the reduced inten- 
sities, we are led to conclude the real space features that 
give rise to this scattering do not carry the periodicity 
of the underlying crystal structure. Four fold rotation 
symmetry is however observed. These facts suggest the 
involvement of an aperiodic interstitial site. The inter- 
stitial Fei site (Fig. (Hb)) is located at the center of the 
primitive square Fe planar unit cell at roughly the same 
distance d = zc from the Fe plane as the Se/Te atoms 
(zw 0.30(3)) [li,[ll. 

Because of the so-called phase problem and to take 
into account other knowledge of the chemical structure, 
we use least squares fitting rather than a direct Fourier 
transform to obtain the real space spin configuration 
from the diffuse scattering. The parameters are mag- 
netic dipole moments for the interstitial site and a to- 
tal of 11 non-equivalent surrounding sites. The num- 
ber of free parameters is reduced by forcing identical 
spin configurations in the two planes sandwiching the 
interstitial allowing them to differ only by an attenu- 
ation factor, 77, to account for weaker coupling to the 
more distant plane. Ordered by distance from the in- 
terstitial site, the distinct dipole moments are denoted 
TO„, where n = 1,2,3, • • -,11. The corresponding dis- 
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placement vectors from the interstitial site are labeled 
r^j, where j indexes symmetry related sites. The pa- 
rameters are inferred by minimizing the least squared 
deviation between the corresponding scattering function: 
5(Q) cx |mo + Eni™"exp(iQ-r,y)(l-|-77exp(iQ-c))|2 
and the observed wave- vector dependent T— difference in- 
tensity in the (HKO) and (HOL) planes. Here toq is the 
interstitial dipole moment. 

The best fit S{Q) is shown in Fig. Wic)-{d). That we 
are able to reproduce diffraction throughout the (HKO) 
and (HOL) planes with a value of z = 0.23(6) consistent 
with structural data, and 77 — —0.16(9) indicating anti- 
ferromagnetic correlation between adjacent planes, con- 
firms interstitial magnetism. The inferred spin configu- 
ration is depicted in Fig.UJ^a). The interstitial dipole mo- 
ment is indicated by the central yellow dot and the sur- 
rounding moments are represented by yellow/blue dots 
(parallel/anti-parallel with mo) - their magnitude pro- 
portional to the area of the dots. We see the near- 
est neighbor (NN) moments are parallel to the intersti- 
tial Fe moment. This is consistent with the acute Fe- 
Se/Te-Fe/ bond angle (60.5°, Fig. [Hb)), which is ex- 
pected to yield ferromagnetic (FM) superexchange [l6| . 
Next nearest neighbors (NNN) on the other hand are an- 
tiparallel to the interstitial moment as expected for the 
obtuse (123.9°) Fe-Se/Te-Fc/ bond angle. Comparison 
to atomic displacement discovered through diffuse x-ray 
scattering from Fci+yTe [l3 shows FM (AFM) corre- 
lated spins are repelled (attracted), which is consistent 
with magneto-elastic displacements that enhance mag- 
netic exchange interactions. FM square plaquettes seen 
along the diagonal direction in Fig. HJa) are a resilient 
feature of magnetism in the 11 series that has also been 
noted in Fei iTe [ll| and Fe- vacancy ordered Kj,Fe2-a;Se2 

We also adjusted an overall spin space anisotropy pa- 
rameter resulting in a ratio of 0.81(13) between 
and S^^. Consistency with the polarized beam value 
of 0.67(12) affirms the elastic T— difference scattering is 
magnetic. Absolute normalization of the intensity data 
further allows extracting ymg — 0.22(3) ^ig. For compar- 
ison the product of the nominal and EDX determined in- 
terstitial density and the squared free ion dipole moment 
of Fe"^"'', mo = 5 /is, consistently yields yrriQ = 0.25 /i^. 

The interstitial together with the two nearest neigh- 
bors is sufficient to reproduce the major features observed 
in the (HKO) and the (HOL) scattering planes. The finer 
details of Fig. [Sfc) however, are obtained only when mo- 
ments beyond reach of direct superexchange interactions 
are included in 5(Q). These display an oscillatory be- 
havior reminiscent of a Friedel oscillation. A naive fit 
of the magnetization density to the r— dependent RKKY 
exchange function yields a characteristic wavevector that 
matches the magnitude of the (^, ^) nesting vector. 

For a more rigorous analysis that links the oscillatory 
magnetism to the Fermi surface structure of itinerant 
electrons as for the charge density in Friedel oscillations, 
we use a five band model with exchange interactions to 
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FIG. 4. Magnetic cluster surrounding interstitial Fe; (a) in- 
ferred from the observed diffuse scattering pattern in Fig. [5] 
and (b) calculated from a five band theoretical model. Yellow 
(blue) moments are parallel (antiparallel) to the interstitial 
and the dot areas are proportional to the moment sizes. 



the two nearest Fe spins. The Hamiltonian consists of 
three terms: 



1~L — T~Lq -\- l~Li, 



"Hi, 



(1) 



where "Hg describes the band structure within the five d- 
orbital model and Hint includes the intra- (inter-) orbital 
repulsion U (U'), Hund coupling J^, and inter-orbital 
pair-hopping, G2 [2S[2l|: 



'Hi, 



u' 



Jh s 



Here i and p. are site and orbital indices, respectively. 
We use a primitive unit cell containing one Fe site with 
the Brillouin zone indicated in Fig. [5] Wave vectors in 
this unfolded zone are denoted by dimensionless vectors 
k — Cla/\/2. We index k = /c^jX -I- ky-y in a coordinate 
system rotated by 45° compared to that used for Q = 
Ha* + Kh* so that = {H + K)t: and ky = {H - K)n. 
The bare static susceptibility is 



Xup,^A(q) = 



dk 



{2n 



Gf^^ik + q, ujn)Gxp{'k, -iu„) 



where a;„ — (2n + 1)ttT and G^^(h,uj) is the orbital 
Green's function. The non-zero elements of the 5^ x 5^ 



interaction matrix are denoted Va 
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V^/ij^y = Jh, and V^v^tj. = G, where ^ v. Within 
the random phase approximation, the full spin suscep- 
tibility is x^(q) = hY.i.vX^^^vi'^. where x^^^ = 
— To describe Fe-planes near a magnetic 

instability with a minimum of parameters we choose 
U = 0.95 eV, Jh = G2 = O.OSJJ, and U' = U - 2Jh- 

T~i\nip describes the exchange interaction between the 
interstitial Fe and neighboring Fe sites: 
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ieNN 



J'S 



jeNNN 



(2) 



Here J < ( J' > 0) is the FM (AFM) exchange constant 
between the impurity spin and the four NN (eight NNN) 
spins in the Fe plane. Himp is treated as a perturbation 
to Hint, with the impurity spin fixed. To leading order, 
we obtain 



s(k) = -4x^(k) 



JCOs(y)cOs(^) 



(3) 



+ J ( COs(y)cOs(^) +C0S( — )C0S(^ 



The structure factor, including the contribution of the 
impurity spin is S{Q) oc |1 + s(k)p. While has nest- 
ing peaks at k = (7r,0) (Q = (i, i)) 



!|, these are 

suppressed by the square bracket in Eq. [3] The fit to the 
experimental data gives J = — 70 meV and J' = 40 meV. 
Consistent with the effective nature of Himp, there is a 
considerable robustness to the fit: the essential features 
are the FM J versus AFM J' and 0.2|J| < J' < 0.8|J|. 
For comparison, the dominant NN and NNN exchange 
constants in Fei.osTe, with similar Fe-Te-Fe bond angles 
are J = - 51(3) meV and J' = 22(4) meV 

The calculated structure factor, iS(Q), is shown in 
Fig.[He)-(f) and corresponding real space magnetization 
map in Fig. |4l[b) . Comparing to the experimental data 
(Fig.lU^a)), there is reasonable agreement up to the third 
NN beyond which the theory overestimates the magni- 
tude of induced magnetization and modulation along the 
(^, ^) direction. Possible reasons include lack of orbital 
specificity to the interaction parameters and effects from 
neighboring interstitial sites. Indeed, the appearance of 
nominally elastic diffuse magnetic scattering in our ex- 
periment indicates a spin-glass like state that links inter- 
stitials. Further information about associated spin dy- 
namics was recently provided for Fe1.01Teo.72Seo.28 |23| . 

Confirming indications from resistivity measurements 
[m and predictions from Density Functional Theory [^J] , 
our data show the interstitial site develops a full lo- 
cal moment. Superexcliange interactions further enforce 
FM plaquettes around impurities with fairly large mag- 
netic moments. Sprinkled at random through the sample, 
these favors spin configurations where the primitive unit 
cell carries dipole moment so that the (^0) type dou- 
ble stripe structure emerges as a compromise between 
the (5^) semi metallic nesting instability and FM su- 
perexchange interactions. Indeed this manifests in our 
impurity band structure calculation (Figj^Je)). 
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FIG. 5. (a) T-dependence of neutron scattering intensity at 
Q=(i,0,0) and huj = 0.0 meV (red), and Q=(i,i,0) and 
hoj = 0.5 meV (blue), (b) DC susceptibility measurement at 
30 Oe (blue diamonds) showing diamagnetic screening which 
yield an upper bound of 92.9(7)% on the superconducting 
volume fraction. Specific heat data (red circles) from which 
a volume fraction of 83% is extracted. 



We now examine the interplay between interstitial 
glassy magnetism and superconductivity. Fig.[5ja) shows 
the T— dependence of inelastic scattering at Q = (^, ^) is 
precipitously suppressed for T < Tn as the s± gap opens 
and the spin resonance develops "^l. The elastic scatter- 
ing at (^,0) on the other hand grows upon cooling with 
no apparent anomaly at Tc- Despite the % level intersti- 
tial concentration, the spatial extent of the associated 
Friedel oscillation (Fig. 2]) ensures microscopic coexis- 
tence with the > 80% superconducting volume fraction. 
The large energy scales (—J, J' >> ksTc) that control 
the interstitial polaron and the different characteristic 
wave vectors associated with magnetism and supercon- 
ductivity are surely relevant here. At the same time pre- 
vious studies show interstitial iron does reduce the super- 
conducting volume fraction [Tlj . These facts suggest two 
length scales are involved as in the mixed phase of a type 
II superconductor: A normal polaron core with a Friedel 
oscillation extending into the superconducting bulk. 
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